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We present simple, concrete, two-fermion models that exhibit thermodynamically stable isotropic 
translationally-invariant gap less superfluid states (breached pair superfluidity). The mass ratio be- 
tween the components and the momentum structure of the interaction are crucial for determining the 
stability of such states: Idealized, momentum-independent ("contact") interactions are insufficient. 

PACS numbers: 64.10.+h,65.90.+i,71.27.+a,74.25.Dw 
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INTRODUCTION 

Recently there has been interest in superfluid fermion 
systems where there exist superfluid states that retain 
gapless fermionic excitations flSSIJl- These states em- 
body "phase separation in momentum space" : some de- 
grees of freedom pair, forming a superfluid, while others 
remain unpaired, maintaining the properties of a Fermi 
surface. They are likely to become experimentally acces- 
sible in the near future []|. 

The important new result is that such states are sta- 
ble, but to see this, one must perform a careful analysis of 
the momentum dependence of the interaction. Appropri- 
ate systems (such as two-band models 0,0) have been 
studied, but these states were missed because the pos- 
sibility of momentum dependence was ignored. We also 
address the questions of instability toward phase separa- 
tion || and local currents [|J raised in response to the 
proposal Q. 

We shall consider systems of two species that, in the 
absence of interactions, would have two distinct Fermi 
surfaces. Simple heuristic considerations suggest the pos- 
sibility that pairing takes place about the Fermi sur- 
faces, but that there is no pairing in a region between 
the surfaces: this led to the term "breached pair" 
A breached pair superfluid state (BP) is characterized 
by the coexistence of a superfluid and a normal com- 
ponent in a translationally invariant and isotropic state. 
These components are accommodated in different regions 
of momentum space with the normal component residing 
in the "breaches" , bounded by gapless Fermi surfaces. 

A state of this type was also considered by Sarma • 
He considered the case of a superconductor in an external 
magnetic field, and found that, although there is a self- 
consistent mean-field solution with gapless modes, it is 
unfavored energetically to the fully gapped BCS solution. 
Similar results were considered in the context of color 
superconductivity 01 > again concluding that these states 
are not stable at fixed chemical potential. 

Since the fully gapped BCS solution enforces equal 
numbers of each species, it was incorrectly suggested Q 
that one might stabilize the gapless phase by enforcing 
constraints on the particle numbers. Indeed, by enforcing 



unequal numbers of each species, one forbids the forma- 
tion of a fully gapped BCS state, but admits "breached 
pair" states in which the excess in one species can be ac- 
commodated by the breach. In the QCD context, a sim- 
ilar argument has been made by imposing charge neu- 
trality [31 ML Recently, however, Bedaque, Caldas and 
Rupak |g, |l2| pointed out that a spatially mixed phase 
may be energetically preferable: this rules out the first 
possibility 2] but may not affect the QCD case due to 
the long-range gauge interactions. 

Here we clarify, broaden, and correct this discussion. 
We conclude that: 

• For extensive systems, one can not stabilize a state 
by imposing different global constraints (such as 
fixed particle number): The composition of the 
state can be completely determined from an anal- 
ysis of the grand canonical ensemble. The specific 
examples considered in |2J are accordingly unstable. 

• With the proper momentum structure, however, 
one may realize breached pair superfluidity in states 
that are thermodynamically stable for fixed chem- 
ical potentials. We exhibit these below. 

Our considerations do not apply directly to non-extensive 
systems. Charge conservation or color neutrality con- 
straints enforced by long-range gauge forces might stabi- 
lize BP phases. (Of course, the possibility of a competing 
mixed phase must still be considered quantitatively.) 



THERMODYNAMIC STABILITY 

In the context of two component fermionic systems as 
considered in 0, E3 > three competing homogeneous 
phases have been considered: a normal state of free 
fermions (N), a fully gapped superfluid phase (BCS), and 
a gapless BP phase. The BCS phase has complete pairing 
between the two species, and thus enforces equal densi- 
ties. The other phases admit differing densities. 

Upon solving the self-consistency conditions (gap equa- 
tions), one commonly finds that over a range of chemi- 
cal potentials there are three distinct solutions. To de- 
termine which is stable in this grand canonical ensem- 
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ble, one must minimize the grand thermodynamic poten- 
tial (equivalently, maximize the pressure) of the system. 
Typically, two of the three solutions are minima on ei- 
ther side of the third BP state which is a local maximum: 
Fig.Oshows a typical potential. The gapless states found 
in |2j, |l0( correspond to local maxima, thus the compet- 
ing state with larger gap parameter A has higher pressure 
and renders the BP state unstable in this ensemble. 

If the stable solution is fully gapped, then it has equal 
densities. By fixing the particle densities to be unequal, 
one may forbid this BCS state. Furthermore, upon com- 
paring the Helmholtz free energies H — which must be 
minimized in this ensemble — one may find that the "un- 
stable" BP state is favored over the normal state N. 

This apparent contradiction in the stability analysis 
based on different ensembles can be resolved by con- 
sidering a mixed phase 0, Il2l | which has an even lower 
Helmholtz free energy H . That such a resolution is al- 
ways possible, however, may not be apparent; indeed, it 
is generally hard to determine the mixed phase explicitly. 
By using general properties of the grand thermody- 
namic potential Q, however, one can argue that such a 
solution is always possible, as follows. By definition, 



PV = Sl(jt) = mm(H — ft ■ N), 



(1) 



minimized over all competing phases. Thus, f2 is a con- 
cave function of the chemical potentials ft = (fi a ,fif,)- 
(We consider here fixed T = 0, but concavity in T also 
follows from maximizing entropy.) Furthermore, there 
is a one-to-one correspondence between tangents Q and 
states of fixed particle number: 
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(2) 



When Vi is not differentiable, there is a cone of possible 
tangent hyperplanes which contact f2 and which bound 
fi from above (see Fig. |T). This cone of tangents de- 
scribes various possible mixed phases composed of the 
pure phases (where £1 is differentiable) that intersect at 
the singularity. To find the state that minimizes H for 
some fixed constraint N = N one simply forms the hy- 
perplane with gradient Nq and drops this until it contacts 
the surface f2. The first point of contact will define ei- 
ther a pure or mixed state which satisfies the appropriate 
constraints. Since this state also lies on O, it minimizes 

for the fixed chemical potentials defined by the con- 
tact point. No matter what constraints we apply, there 
is always a stable state in the grand canonical ensemble. 

This argument is valid only for extensive thermody- 
namic systems. Long-range interactions can render the 
energy of some pure phases non-extensive (due, for ex- 
ample, to the rapidly diverging Coulomb energy per unit 
volume V as V — > 00). In such cases, a mixed phase 
would contain bubbles of limited size. The surface en- 
ergy of these phase boundaries becomes a volume effect 



and must therefore be taken into account, even in the 
thermodynamic limit (see for example |l4|l. This com- 
plicates the relation between N and fx. In this rest of this 
letter, we shall consider only finite-range interactions. 




FIG. 1: The cone of tangent (hyper)planes to a thermody- 
namic potential density — P = Q.(p)/V. Immediately to the 
left of [io is a pure phase with density ul while immediately 
to the right is another pure phase with density ur. The den- 
sities are the negative slopes of the tangents at fio according 
to ©. At /1 = fio there is a continuum of mixed phases: 
These consist of a volume fraction x at density ul and the 
remaining fraction 1 — x at density ur. The average density 
over all space, n = xn^ + (1 — x)nu, lies within n € (til, tir). 



STABLE BREACHED PAIR SUPERFLUIDS 

We now demonstrate, by example, how to realize pure 
BP superfluid states in extensive systems. We shall con- 
sider the mean-field analysis of two models, each with two 
species of fermions a and b of differing masses m a < m&: 
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Hi. (3) 



The models have different interactions TLf. and 
We shall consider these systems in the grand canonical 
ensemble at zero temperature by minimizing the thermo- 
dynamic potential density fi(/i a ,//(,)/V\ It will be nat- 
ural, however, to use the parameters pf = ^2m,i^i in 
place of the chemical potentials 

The first model posits a spherically symmetric static 
two-body potential interaction V(r) between the two 
species a and b: 
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d 3 fd 3 y V(\x-y\) atblb ff a s . 



(4) 



Defining m± = 2m a mb/(mb ± m a ), Li± = (fj, a ± A*b)/ 2 , 
and 
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and considering only homogeneous (translationally in- 
variant) and isotropic phases, we find that extrema of 
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FIG. 2: Qualitative T = phase diagram for interaction gj 
with a Gaussian potential V(r) oc exp(— r 2 /2A 2 ). All mo- 
mentum scales are in units of ft/ A and all energy scales are in 
units of ft 2 /(m+A 2 ). The mass ratio is mb/m a = 50 and the 
coupling strength has been chosen so that 2m+A PO /po = 0.1 
at the point marked "x" where (pi, ,p a ) = (13,10) to en- 
sure weak-coupling. (This ratio is less that 1 throughout this 
diagram.) Note that the lower BP region has more heavy par- 
ticles b while the upper BP region has fewer heavy particles. 
The upper type may be realized in the QCD context Q- All 
phase transitions are first order. 



satisfy the gap equation 
d 3 q 
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2 J (4)* + A2 
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where V(p) is the Fourier transform of V(r). The in- 
tegral © runs over the region R outside any "breach" . 
R contains momenta where the two quasiparticle disper- 



sions Ep 



(7) 



have opposite sign. (See |2( for further details about the 
generic breach structure.) Note from (JBJ that A p is gen- 
erally largest about po, where e+ = 0. 

Equation © can be solved numerically to find ex- 
tremal points of the thermodynamic potential. Over this 
set of self-consistent solutions, one can minimize Q to 
determine the phase structure. 

We have done this for a variety of interactions, and 
find similar qualitative structure: a central strip of fully 
gapped BCS-like phase about p% — p£ , with normal un- 
paired phases outside (see Fig. [SJ) Depending on the 
model parameters, these phases may be separated by a 
region of BP superfluid phase. To verify that these in- 
deed contain gapless modes we plot in Fig. [31a sample set 
of occupation numbers, quasiparticle dispersions, and the 
gap parameter A p , 
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d 3 g- 
(2V)3 



V(\p-q\)(bpa-p). 



(8) 
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FIG. 3: Quasi-particle dispersions (top), occupation 
numbers n a and ni (middle), and gap parameter A p for a 
sample BP state at (pf,pf) — (13,10). All units and pa- 
rameters are described in Fig. |2] Notice that there are two 
"Fermi" surfaces at p » 10.3 and p = 13. The first occurs 
where A p becomes too small to support the gap. The second 
is simply the Fermi surface for b which is virtually unaffected 
by the pairing. The "breach" occurs between these surfaces 
and only the region R outside contributes to the gap equa- 
tion ©. 



The presence of gapless fermion modes depends cru- 
cially on two factors: 1) the momentum structure of A p 
and, 2) the mass ratio. 

First, A p must be large in some regions and small in 
others. If A p is large enough at a Fermi surface, it will 
induce pairing at that surface and support a superfluid. 
If it is also small enough at the other Fermi surface, it 
will not appreciably affect the normal free-fermion be- 
haviour. The problem with previous analyses is that 
they assume pointlike interactions, implying a constant 
A p = A. Physical interactions, however, tend to exhibit 
more complicated behaviour and suitable A p are quite 
generic: A p tends to peak about the Fermi surface of the 
lighter species and fall off to at least one side. The model 
shown in Fig. |3 and Fig. |3 for example, has a Gaus- 
sian interaction. Longer-range forces (such as a screened 
Coulomb interaction) tend to plateau to the left of po but 
still fall to the right of pq. 

Second, as was emphasized in the one may reduce 
the cost associated with shifting the Fermi sea pF by in- 
creasing the mass m&. Thus, by choosing a large enough 
mass ratio, one may always move the Fermi surface for 
the heavy species to region where A f is small enough to 
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FIG. 4: T = phase diagram for model © with a hard 
cutoff f(p) ~ 6(p — 10a) that has been smoothed over the 
range from p £ (9.7, 10.3). All momenta are expressed in 
units of a where 10er is the cutoff scale, and all energies are 
expressed in units of <7 2 /(2m+). The mass ratio is mb/m a = 4 
and the coupling g has been chosen so that 2m+A/po = 0.2 
a t Pa — Pt — Po — 10cr to ensure weak-coupling. (This ratio 
is less that 1 throughout this diagram). All phase transitions 
are first order as discussed in the text. The mixed phases 
of sua would be found on the solid lines. The sample state 
in Fig. at (pf,pf) = (11.5,9.2) is marked "x". 

leave the Fermi surface undisturbed. The states shown 
in Fig. and Fig. [3] have a mass ratio mb/m a = 50. 

Since the variational states of model (@J are parame- 
terized by a variable function A p , the set of states over 
which the minimization Q must consider is enormous, 
and we cannot be certain to have found the global min- 
imum. We have searched for stable fixed-points of the 
gap equation JfjJ) and compared them, so our results for 
this model are consistent and plausible, but not rigorous. 

The second model allows us to be rigorous. We fall 
back to the type of factorized, cutoff interaction often 
considered in BCS models: 

Ul=9 J ff3 ( ffs /(P)/(g)^^ S -f^ ^ 

In this model, the combination Af(p) plays the same 
role as A p in first model, with the cutoff function f{p) 
providing the required momentum dependence. The trial 
states are now parameterized by a single number A = 
(A) where A = g J d 3 p/(2ir) 3 f(p)b p a_ p . One can now 
find the global minimum by plotting (see Fig. |5J) 

= min (0\n - fi a h a - fi b n b \0) , (10) 
V (e\A\e)=A 

where we minimization over all BCS style ansatz \9) with 
given expectation A. This minimization is equivalent to 
comparing all solutions of the mean-field gap equation 




We conclude that, within the mean- field approximation 
of homogeneous phases at zero temperature, this model 
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FIG. 5: Quasi-particle dispersions Ep (top), and occupation 
numbers n a and nj, (middle), in a sample BP state. This state 
has gap parameter A w 11 which is the global minimum of the 
grand thermodynamic potential density Q(A) /V (bottom) as 
defined in JTOJ. The maximum at A w 5.6 corresponds to 
an unstable BP state. These figures correspond to the point 
(p[,Pa) = (11.5,9.2) in Fig.H 

has the phase diagram shown in Fig. 21 We plot the 
properties of a sample BP state in Fig. \5\ to illustrate 
that there are indeed gapless modes. 

To model A p more accurately one might use a func- 
tion f(p) where the location of the cutoff stays near 
Po. This introduces an inconsistency in the thermody- 
namics because f(p) is really a property of the Hamil- 
tonian, while po depends on the chemical potentials fa, 
thus N ^ —dCl/dfj,. For small coupling and high densi- 
ties, these spurious dependencies become small and the 
resulting phase diagram is qualitatively like Fig. [21 

Finally, we address the issue of the instability discussed 
in |jj where they claim that the superfluid density is 
negative due to a large negative contribution from the 
diverging density of states at E = near the transition 
to the BP state. If one simply computes d 2 f2/dA 2 , one 
finds exactly the same negative contribution indicating 
that the BP solution under consideration is an unstable 
maximum rather than a stable minimum. The solutions 
we present here are all global minima, and hence stable. 

This raises an interesting point: if the BP/BCS transi- 
tion were second order, then the density of states would 
formally diverge. Indeed, one finds d 2 il(A)/dA 2 = — oo 
at certain chemical potentials. Near the BP/BCS tran- 
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sitions, fi(A) develops a cusp separating two competing 
local minima: one is BCS and the other is BP. Thus, 
in the T — mean-field approximation, the transition 
must be first order. At finite temperature, the cusp is 
smoothed and we suspect that the transition line ends 
at a critical point. In this way, the T = BP transition 
avoids instability. In non-extensive systems such as QCD 
where gapless states may be stabilized by neutrality con- 
straints, similar instabilities have been noted The 
resolution may be the formation of a non-homogeneous 
phase. This possibility requires further analysis. 

Realizations of a stable BP phase require either non- 
extensivity, or a finite-range momentum dependent in- 
teraction with a large mass ratio. The former may occur 
in high-density QCD 0, 0] where gauge interactions may 
stabilize the state. The latter may occur in a quantum 
gas of cold neutral atoms operating near Feshbach reso- 
nance with effective masses tuned by a laser lattice M, in 
a system of trapped ions with dipolar interactions 
or in superconductors with overlapping bands 0, @ ■ 
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